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Abstract 

The renormahzation method of Bogoljubov-Parasiuk-Hepp-Zimmermann 
(BPHZ) is used in order to derive the renormahzed energy shift due to the 
%^J . gauge invariant Kallen-Sabry diagram of the two-photon vacuum polariza- 

If^ \ tion (VPVP) as well as the self energy vacuum polarization S(VP)E beyond 

the Uehling approximation. It is outlined, that no outer renormalization is 
required for the two~photon vacuum polarization and that only the inner 
f^ I renormalization has to be accomplished. It is shown that the so-called non- 

gauge invariant spurious term is absent for a wide class of vacuum polarization 
(VP) diagrams if one applies the widely used spherical expansion of bound and 
free-electron propagator. This simplifies significantly calculations in bound 
state quantum electrodynamics. As one result of our paper the use of the 
BPHZ-approach in bound state QED is established. 
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, J 1 Introduction 

X 

d ' Highly charged ions provide an ideal scenario to demonstrate the validity of QED 

in strong electric and magnetic fields by measurements of the Lamb shift at utmost 
precision. In this respect the recent experimental progress made in measurements 
of the Lamb shift in hydrogenlike ions |jl|, indicate that calculations of all QED 
corrections of order a^ {a ~ 1/137.036 is the fine structure constant), but exact in 
the coupling constant Za become relevant. The present status of theoretical predic- 
tion for the Lamb shift have been presented in [0, |^ . Even the most difficult gauge 
invariant set, the second-order electron self energy correction, have been recently 
calculated for the hydrogenlike ions uranium and bismuth |^. The renormalization 
of these diagrams has been carried out in Ref. [^ and within a different approach 
in Ref. 0| . At present two of these two-photon diagrams are known in the Uehling 
approximation only, namely the two-photon vacuum polarization VPVP and the 
effective self energy S(VP)E (see Fig. |T|). Concequently, mainly the VPVP cor- 
rection remains an additional source of theoretical uncertainty in the Lamb shift 
predictions of hydrogenlike ions. P, H] . The main contribution originates from the 



Uehling part of the VPVP diagram (—0.6 ± 0.12 eV for uranium and —0.34 ± 0.07 
eV for lead). To provide a result, complete in order o? ^ for the Lamb shift is in- 
deed a challenge for theory. As a step towards this goal we derive the renormalized 
expression for the two-photon vacuum polarization. At the same time we derive a 
renormalized expression for the higher-order energy shift of the effective self energy 
diagram S(VP)E. This diagram is of special interest not only since its contribution 
to the Lamb shift is unknown for all elements in nature but also due to the fact 
that it contains the one-loop vacuum polarization (VP) as a subdiagram. It is well- 
known that in general the evaluation of VP is connected with a so-called nongauge 
invariant spurious term. We shall show that this spurious term is absent for a wide 
class of diagrams which contain the VP diagram. This result simplifies essentially 
evaluations of such diagrams in bound state QED. 

The paper is organized as follows. In section 2 we explain the BPHZ-renormalization 
approach and it's application to bound state QED. The renormalization of the VPVP 
diagram is discussed in section 3. It is shown that no outer renormalization is re- 
quired and the renormalized expression is specified in section 4. In section 5 we 
shall discuss some basic properties of the fourth-rank vacuum polarization tensor. 
In section 6 the renormalized expression of higher-order of the S(VP)E diagram is 
derived. Finally it will be demonstrated that the nongauge invariant term is ab- 
sent if one applies the partial wave decomposition of the electron propagator. The 
renormalized expression of S(VP)E diagram is specified in section 7. 

2 BPHZ— Renormalization scheme for QED 

The general renormalization scheme of QED is formulated for graphs with free- 
electron propagators. Therefore, in order to isolate the ultraviolet divergencies 
we first have to perform a potential expansion [|I0], |Tl], |12| of the bound electron- 
propagator 

+ I rfV I St' SUri - r, cohoVir)Spir, r', cu^oV (r') SUr' - r^, cu) (1) 

that is represented graphically in Fig. 0. After this it is possible to apply the stan- 
dard renormalization prescription for the free QED. We use the BPHZ-renormalization- 



method |T^, |14[. For any given divergent loop-integral Fr corresponding to a Feyn- 
man diagram F this approach allows for a systematic isolation of the divergent parts 
and to derive the finite contribution Fp of the loop-integral under consideration. 
This procedure is formally expressed in terms of Bogoljubov's 7^-operation acting 
onto the integrand Jp of the divergent loop-integral F^. This operation determines 



the integrand Rr, which yields the finite integral Fp: 

Fr= f d% ... d% /r -^ F; = / d\ ... d% Rr . (2) 



Depending on the superficial degree of divergence ci;(r) of the diagram F the action 
of Bogoljubov's 7?.-operation is defined recursively via 

^^^iii-f)% .f .(r)>ol 

[ Rr if cu(r) < J 

where Bogoliubov's recursion formula is given by 

k 

% = Ir+ E ^m.u...,,.}Ili-t'nK ■ (4) 

Here ^(7) = A — Eb — \Ep defines the superficial degree of divergence for QED where 
Eb and Ep denote the number of external bosonic and fermionic legs, respectively, 
asigned to an arbitrary diagram or subdiagram 7. The sum in eq. (^) runs over all 
combinations of disjunct superficially divergent sub diagrams 7^ of the diagram F, 
and 7^ denote the superficially divergent subdiagrams of the subdiagrams 7j. The 
mathematical expressions for the subdiagrams are obtained by means of the usual 
Feynman rules. The symbol Ir/-yi,...,-yk stands for the integrand of the diagram F 
where the subdiagrams 71,.. .,7^ are merged to a point. The operator V denotes 
the Taylor expansion with respect to the independent external electron momenta 
Pi and photon momenta kj up to order cc;(7) of the diagram 7. According to the 
renormalization group equation of QED the reference point of this Taylor expansion 
can be choosen arbitrary. However, a special choice of a renormalization scheme 
with fixed reference point implies renormalization conditions which are kept fixed. 
Here we employ on-shell-renormalization conditions and utilize the gauge invariance 
of the coupling between electron-positron field and the radiation field. This implies 
the following conditions for the renormalized electron self energy, the renormalized 
vertex function and for the second- and fourth-rank vacuum polarization tensors of 
the photon, respectively, in perturbation theory of order (a)" (n = 1, 2, 3, ...): 



P=m 






nW'-^-W =0, fiW;r(fci,fc2,A:3 



iiupa 
k=0 



= , (5) 

ki=k2=k3=0 

where m denotes the physical electron mass and ki = means actually ki = 
(0,0,0,0). Therefore from the renormalization-conditions (^ one can define the 
Taylor-Operator P up to first order according to the prescription: 
The operator f^ expands the integrand belonging to the diagram 7 



a) with respect to all independent Np external electron momenta ^^ around m 

b) with respect to all independent Nb external photon momenta kj around 0. 

In our context here it will be sufficient to perform the Taylor expansion up to first 
order in external momenta. In mathematical terms we can write 



p=m,k,=0j 




+ ... . 



3 Renormalization of VPVP 



(6) 



The contribution of the Kallen-Sabry-diagrams with the free electron propagator 
are well-known p, 0. Accordingly we subtract these terms from the two-photon 
vacuum polarization diagram and concentrate on the higher-order contribution of 
VPVP which is graphically shown in the Fig. |^. As the major result in this chapter 
we derive the corresponding renormalized energy shift. 

3.1 The outer renormalization 

According to the BPHZ-renormalization scheme as defined above the case where 
the Taylor-Operator t^ acts on the whole diagram F may be called outer renor- 
malization. From Eq. (^ it becomes evident that outer renormalization is always 
required for a superficial degree of divergence uj{r) > 0. In the case of the VPVP 
diagram we have the following situation: 

In the potential expansion of the energy shift in Fig. ^ diagrams with at least four 
external Coulomb-photon legs occur. According to the Furry's theorem diagrams 
with an odd number of external photon or bosonic potential legs vanish. The su- 
perficial degree of divergence u is zero for diagrams with four external bosonic legs 
and is negative for all diagrams of higher-order in (Za)" with n>4. In view of Eq. 
d^) outer renormalization is required if and only if the number of outer bosonic legs 
is equal to four. These diagrams have to be considered with special care. 
At this point we note that in the case of the one-loop vacuum polarization VP one 
must consider the diagram with four outer bosonic legs with special care as well 
(see for example [|16|]). Employing any regularization scheme in the one-loop case 



VP the corresponding counter term vanishes. Moreover, if the spherical expansion 



for the free and bound electron propagator is used no spurious non-gauge invariant 



term occurs in the 1-loop case VP IT^- For detailed discussions we refer to the 
literature ([0, |T^, ^ and [0). With this in mind one could expect that, after any 
regularization is performed, the counter terms required for the outer renormaliza- 
tion may vanish in the two-loop-case VPVP as well. In any case an explicit proof 
is desirable. 

We turn to the investigation of the couterterms of the outer renormalization for the 
VPVP diagram. Performing the potential expansion of the energy shift shown in 
Fig. 1^ leads to six diagrams with u = (see Fig. |). To give an explicit example 
we consider graph Pa depicted in Fig. |^. This diagram includes the subdiagram 7a. 
It's contribution to the energy shift of an arbitrary bound electron state ifn reads 

XV(*=0 V(k.) V(k,) o._(fe. + fc^ + ,3p + ,, O'o FrJkuk^M) ^.M . (7) 

where V{kj) {j = 1, 2, 3) is the Fourier-transform of the Coulomb potential of the 
nucleus: 



V{kj) = I d^re~'^^''V{r) . 



The explicit mathematical expression for the Feynman diagram is given by (with 
h = {Q,ki) ,^ = 1,2,3) 

d^qi d^q2 



Fr^{ki, ^2, fcg) = J ,^^|^ .^^,^^ ^rA(?i, 12, ^i, ^2, ^s) , 



(27r)4 (27r) 
IrAiQu (l2, ki, k2, ks) 



xiTi 



1 

7o^ — —Ti — —r, — —r, —^0-, — —Ti — —Ti —lo 



^^ + ^-^^ + It^ + It.^ - m + ie ^^ + It^ + ^t^-m + ie ^^ + ^t^ - m + ie 
1 1 „,\ 1 1 



X -ie"^ . 7c,- — -—7" I / , ,, -— 7o 



gl + ^e ^1 ~ ^2 + 1^3 ~ "^ + ^^ J ^^ + It^ — m + ie ^^ — m + ie 



(9) 



In round brackets the mathematical expression of integrand refering to the subdia- 
gram 7a (one-loop-self energy) has been inserted. Applying the 7?.-operation Eq. 
d) and (ID yields 

Rr^ = ^ji^^^ + ^r°f ^■- (10) 



with 



RrT = ^r. - /r./,. {t'^Q 



RrT = -t'^^r, + 1'^ (^r./,. (t^^4)) 



(11) 



The Taylor-Operator t^^ acts up to zeroth order while the other one P^ acts up 
to first order since uj(Ta) = and ci;(7a) = 1- Accordingly, the terms for the outer 
renormalization are 



and 



(27r)4 (27r)4 



Al A2 Pa 



^FaI^I, 92,^1,^2,^3 



(fqi d^q2 

(27r)4 (27r)4 



/rjgi,g2, 0,0,0) (12) 



(27r)4 (27r) 



^Tf^hMl^ ^1' ^2, ^3) (i^^^^A (9l> ?2, fcl, ^2, ^3 



'^^'?1 f / r> r> r>\ , v-fD' /" ^^'?1 



^^'7(211^^^-^^^'°'°'°^ + S«'/-^(4-m)^^(gi, 52, 0,0,0 



(2vr) 



(13) 



respectively. The denotation (^^ — m) Iva/'ja implies insertion of the Dirac-structure 
(^1 ~ '^) ^'t the same place in the diagram rA/7A where there was the subdiagram 
7a before. In Eq. ([T3|) we used the usual definition of the one-loop-counter terms 
S(^) and S(i)' of free QED: 



s(i) = - 



le 



d% 1 a 1 

(27r)^ °^^ - ^2 ~ "^ + *^ (ji + ^e 



(14) 



and 



^>^'"'--^^/ 



9 /■ d'^q2 



7la 



7 



(271)"^ ^^ - J^2 ^ "^ + ^^ 5'2 + ^e 



(15) 



Similar steps have to be performed in all other diagrams of Fig. ^. With the aid of 
the one-loop Ward-identity 



dq^i 
together with the identity 



dq[ 



'h - ^) ^rA/7A = 7o/rB/7B 



(17) 



and the relations between renormalized and unrenormalized vertex and self-energy 
operator 



S«(^i) = S« + (^1 - m) S«' + S(i 



)rcn / 



111 •> 



all counter terms refering to the outer renormalization can be collected. One ends 
up with the following expression for the sum of all counter terms generated by the 
operator R^^^'^^ corresponding to the six diagrams of Fig. ^: 

u, qi / ti y2 fr) £) outer 1 o £> outer 1 o £> outer 1 £) outer 1 £) outer 1 £) outer\ 



/u, qi / ti y2 fr) £) outer 1 o £> outer 1 o £> outer 1 £) outer 1 £) outer 1 £) 
72^7 72^1^^^- +^^^B +2i?ro +^r, + Rr, +R 



rp y 



Tr 



d qi d 



(27r)4 9g? 
1 



1 t ClVpn . , , < 1 1 



(l)ren 



— -m + «e a-, — m + ze a-, — m + le 



ll 110 I o^ ^^ ,,u I „u 1^^ 



y>(l)ren 



7oi r— 7oi -— 7o 



n 



m + ie ""^ ^-^ — m + ie ^^ — m + it ^^ — m + ie 



(19) 



The last integral is understood as being regularized. The renormalized self energy as 
well as the renormalized vertex correction do not posess any singularities. Therefore 



the total integrand has simple poles only at g° = ^J q\ + {m — ie)"^. Accordingly, 
a Wick-rotation g° — >• iq^ to Euclidean space can be carried out and all surface 
terms vanishes after regularization. This implies that all counter terms generated 
by R-^'^^'^^ cancel each other. Therefore all singularities associated with the outer 
renormalization must cancel in the unrenormalized expression of VPVP. 

3.2 The inner renormalization 

According to the previous subsection only an inner renormalization associated with 
the action of R^'^^^^ is required. For this pupose we perform the potential expan- 
sion of the bound-electron propagator until three outer bosonic potential legs. An 
appropriate potential expansion of the energy correction of higher-order in eq. @ 
is given in Fig. ^. For the last two diagrams C2 and C3 an inner renormalization 
is not adequate. For the other six diagrams Al, A2, A3, Bl, B2 and CI a simple 
inner renormalization is necessary. There are the inner divergencies of the one- 
loop-self energy as well as the one-loop-vertex correction. The renormalization of 
this one-loop-expressions has to be carried out using eq. (JT^). Accordingly, one 
should perform the replacements of Fig. ^ in order to receive the full renormal- 
ized expression for VPVP. If one uses the Ward-identity Eq. (jIBp, Furry's theorem 



and the potential expansion shown in Fig. 0, it is not difficuh to reahze that all 
the terms which are proportional to T,^^^ are cancelled against the terms which are 
proportional to A'^^^ Especially, it holds the following relations for these counter 
terms: 

2A(^) CT(B1) + 2S(^)' CT(Al) = 0, 

A(^) CT(B2) + A(^) CT(Cl) + 2 S^^^' CT(A2) + S^^)' CT(A3) = . (20) 

Finally, if one collects all terms proportional to S*^^^ and retransforms the potential 
expansion one ends up with the renormalized expression for the total energy shift 
of the VPVP, see Fig. |. 

4 The renormalized expression of VPVP 

Now we proceed to write down the mathematical expression of the renormalized 
vacuum polarization correction. The energy shift in higher-order reads: 

^^VPVPrenh.o. ^ J d'v, Tf^in) ^o ^nin) 

(21) 

It is depicted in Fig. ^. The potentials are given by (note that S^^^ is proportional 

to e^) 



oo oo 



U'^^^^ir^) = e* / d-'r2 / dVg / d-'r^ 



27r J 2n 

oo — oo 



xTr [705*^(^1, r2, r3)7„S'F(Ei - E2, rs, r4)-fpSF{Ei, r4, 7-2)] . 

(22) 

dEi f dE2 



U^'ir,) = e'Jd'r2d'r,d'r,d'r, J ^ J 



27r J 27r 

CO —00 



xD'^'ir, - r2, 0)Z^"'^(r3 - 7-4, ^2)1^(^5) 



xTr 



SUr2 - ^3, E.h^SUrs - r,, E, - E2hfsSUr4 - n, E,] 



xSUr6-r2,E, 



(23) 



00 , „ 00 



rBO 



u^'in] 



e* J d'r2 d'rs d'r, d'r, J ^ J 



dEi f dE2 

~2^ 



-00 —CO 



.a/3/ 



xD^'^r^ - V2, 0)Z^"^(r3 - r^, E2)V{r^) 



xTr 



S%{r2 - rs, E^)^^SUr3 - n, E^ - ^2)5^(^5 - n, E, - ^2)7/3 



xSUr4-r2,Ei] 



(24) 



rCTl 



U^^'iri] 



Jd'r2d^r, J _D00(ri-r2,0)Tr [SF{r2,rs,E)joSF{rs,r2,E)] , 



(25) 



f/^'^'^(n) = e' I d'r2d'r,d'r, I ^D''{r,-r2,0)V{r,) 



dE 



xTr 



7o^f(^2 - 7-3, ^)5^(r3 - 7-4, ^)^^(r4 - r2, ^)1 . (26) 



In the expressions above Sp denotes the bound-electron propagator and Sp denotes 
the free-electron propagator, respectively, and D^^'^ is the photon propagator (see 
appendix) . 

5 The fourth— rank vacuum polarization tensor 

In order to derive a renormalized expression for the self energy vacuum polarization 
S(VP)E it is useful to mention some properties of the fourth-rank vacuum polariza- 
tion tensor. Especially we reexamine some features of this tensor which are related 
to the gauge invariance of the vacuum polarization diagram VP. The fourth-rank 



9 



vacuum polarization tensor is defined as amplitude of the light-light scattering which 
is shown in Fig. llOl The expression of the unrenormalized tensor reads: 



t^mTa" (fci, h, h, m) + 1-;- (fci, ^2, -h -h- 2ks, m) 



nrj,;\(fci,^2,fc3,m) 



where the single three tensors T^ 

g- 



(27) 



are corresponding to the expressions belonging 
to the diagrams depicted in Fig. |10| and they are defined by 



TZTA^iM.h.m) 



(2vr)' 



rTr 7m 



1 



1 



1 



-lu 



1 



m + ie ^ + |i^-^ + 1(^2 ~ ""^ + ^^ 



7a 



X- 



7<7- 



+ lf:^ + 1^2 + If:^ — m + ie ^ — m + 



le 



{2i 



For our purposes it is necessary to consider the general case where the external mo- 
menta ki are not on mass shell [kf ^ 0). The general structure of the solution for 
this tensor is complicated and was derived at the first time by Karplus, Neumann 
in |20|. For our intension here it is only important that the single tensors T 



lJ,,u,cr,X 



are 



unren 
fi,u,a,X 



is convergent. This surprising be- 



logarithmically divergent but their sum 11 
haviour of 11^'^^^^ ^ ^^^ been proven a long time ago |21|, ^, see also [^. Therefore 
a renormalization of n""^J^;s^ seems to be not necesarry, which, however, is not true. 
The vacuum polarization tensor 11"°^^;^ is not gauge invariant and therefore a renor- 
malization is necesarry. Following the BPHZ-renormalization method described in 
section 2, we obtain the renormalized expression 

KZ,.,xiki, k2, ks, m) = tlZ:^^,{k„ k2, k„ m) - fl~,(0, 0, 0, m) . (29) 

Gauge invariance implies (see for example p2|, |l^, 0, |2^ and references therein) 

^rn;:::,.,A(^i,fc2,A:3,m) = o, 
fc2n-:i,.,A(fci,^2,fc3,m) = o, 

klW^l^^^^{k,MM,m) = Q. (30) 

In other words the counter term n"'^J'^;^(0, 0, 0, m) ensure the gauge invariance of 
the fourth-rank vacuum polarization tensor. We note that the regularized tensor 

n;i.,A(^i, ^2, fcs, m, M) = ilfZAk,, k2, ks, m) - flZ::,,{ku k^, ks, M) , (31) 

is also gauge invariant (M denotes a large electron mass of the Pauli-Villars-regula- 
rization). A renormalization of the regularized tensor in Eq. ( PT| ) is not necesarry 
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because the counter term n™J'^J^_<^(0, 0, 0, m) would vanish if one performs such a 

regularization. Sometimes the subtracted term nj^^J™;^(fci, ^2, ^3, ^) in Eq. (PT|), 
which ensures the gauge invariance hke the counter term in Eq. (pOf), is called 
spurious term. Similarly we can consider the counter term n"'^J™;^(0, 0, 0, m) in Eq. 
(|29| ) as a spurious term. 
There are regularization methods, in which the counter term n™J^;^(0, 0, 0,?7i) in 

Eq. (H) or the subtracted term tl'^'^l^^^^{ki, fcs, h, M) in Eq. (^ will vanish |l^ |l9l. 
It means that these terms are not necesarry in such cases in order to ensure gauge 
invariance. Because of this it is important to note that both approaches alone, Eq. 
( P^ ) as well as Eq. (PT]), ensure the gauge invariance of the fourth-rank vacuum 
polarization tensor. Therefore is does not play any role which term will vanish. 
However, because of the simplier mathematical structure of the spurious term of 
Eq. (p9D in comparison with the spurious term of Eq. (|31|) it turns out to be 



much easier to use the approach of Eq. (^) instead of Eq. (0). As we will see 
in the next section the spherical expansion of electron propagator is such a special 
condition where the counter term in Eq. (|29D will vanish. 

6 Renormalization of the S(VP)E diagram and 
disappearance of the spurious term 

In order to isolate the ultraviolet divergencies of the S(VP)E diagram we first have 
to perform a potential expansion of the bound electron-propagator according to 
Eq. (|lD as we have seen already in the case of the VPVP. For the diagram of 
the S(VP)E such an appropriate potential expansion is shown in Fig. |TT]. The 
diagrams A3, B2, C2 and D are convergent. The sum of diagrams A2, Bl and CI 
are nothing else but the fourth-rank vacuum polarization tensor. Therefore they are 
separately divergent but their sum is convergent but not gauge invariant as we have 
discussed in the previous section. It should be mention here that only for hydrogen 
an evaluation of the diagrams A2, Bl and CI was performed (see [^ and references 



therein). At present, for all other elements of the periodic system their contribution 
is unknown. Diagram Al is the so-called Uehling contribution of the S(VP)E. This 
Uehling contribution was evaluated for low- and high-Z atomic systems few years 
ago ||2^ and gave the leading contribution at least for high-Z-systems like uranium 
and lead. 

Now we proceed to prove the disappearance of counter terms of diagrams A2, Bl and 
CI if one uses the spherical expansion of the free-electron propagator (see appendix). 
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Their corresponding energy shift is given by: 

r.(A2)+(Bl)+(Cl) 
n 

1 1 

X 



dEi f d?ki f d^ks f d^k^ 



X 



Ef - fci^ + ie El - (fci - fca - ^4)^ + ie En - Ei - E^(l - ie) 



xn^^^^o,o(^i, fci; El, fci + fcg + ^4; 0, fcg; 0, ki] m) . 
Using the free-electron propagator 



5°(E,ri-r2) 



d^k J^ipi{k)e''"'^ipl{k)e-'''''^ 



(277)3^ E-Ekil 



te 



(32) 



(33) 



where the functions <fi{k) are the solutions of the free Dirac equation in momentum 
space: 

{^ - m) ipi{k) = , ipl^{k)ipi,^{k) = 6i,,i 



,«2 ) 



(34) 



(here the indices i indicate the four linear independent solutions of the free Dirac 
equation) one may derive the relation 



y^ y?i(fc2)yl(fc2) 



(35) 



-E'270 - ^27 -m + ie 

The last three equations lead immediately to the following expression for the counter 
term of the fourth-rank vacuum polarization tensor in Eq. (^): 

n^,.,o,o(0,0,0,0,m) 



—00 



xTr 



d^ 



dE^ 



(a/,5'^(E2, ri - r2)a^,Sp{E2, r2 - r^ 



■ (36) 



Using the identity [H^"'^'^ is the free Dirac operator without Coulomb potential) 

a . = (^0°^-^ - e) r, - r,- {h^'^-^^ - e) , j = 1, 2, 3 (37) 

and the relation 

(^o°---E)5°(E,ri-r2) = (5(ri-r-2), (38) 
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and their conjugate version, respectively, it is not difficult to see that the spatial 
components in Eq. (|36D do not contribute. Finally, with ao = 1 we get 



n^=o,i.=o,o,o(0,0,0,0,m) 



1 1 dE, 



6 



27r 



Tr 



dEl 



{SUE'. 



2,ri -r-i 



(39) 



If one inserts the spherical expansion of the free-electron propagator (see appendix) 
one ends up with 



fiti=o,u=o,o,o{^, 0, 0, 0, m) = - ^ 



dEo d^ 



27r J 27r dEl 

— oo 



Gl^,iE,,n,n) + Gf,^,iE2,n,n 



(40) 



The same expression was derived and investigated already in WTl where it has been 
shown that this expression vanishes if the summation over k is restricted to a fi- 
nite number of terms. This can be seen by inserting the explicit expression of the 
spherical expansion where one gets after a Wick-rotation the following expression: 



n;,=o,i.=o,o,o(0,0,0,0,m) 



1 ^ kl 92 



1 



E 



. 00 

\k\ f U,-C/2 



dE. d^ 



27r J 27r dE^ 

—00 



Re U 



GlU^E,) + Gl%{zE,)]) 



6 V 27r dEl 



X 



{zE2 + 1) VI + El j|.+i/2|-i/2(Vl + El n) h\'>_,Jz^l + El n) 



+{^E2 - 1) VI + El j|.-i/2|-i/2(Vi + El n) /^£U/2|_i/2(Vi + ^2 n] 



E2=oo 



E2=-co 



(41) 



Here j\K±i/2\-i/2iz) and ^L±i/2|-i/2(-^) ^^^ ^^^ spherical Bessel function and spherical 
Hankel function of first kind, respectively. In the last line we used the Gauss law 
in Euclidean space. It can readily be seen that the last expression vanishes if one 
performs both derivations. So it has been proven that the counter term is zero if 
one uses the spherical expansion of the electron propagator where one has to sum 
over K in the last step. This result simplifies essentially the evaluation of arbitrary 
diagrams of bound-state-QED which include the one-loop-vacuum polarization. 
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7 The renormalized expression of S(VP)E 

As we have seen in the previous section the counter term vanishes and therefore it is 
only necessary to subtract the known Uehhng contribution of the S(VP)E diagram. 
The renormahzed expression of the energy correction of higher order is shown in 



Fig. |I2| and the corresponding energy shift reads: 

dEi f dE2 f d?ri r (fr2 f d^r^ f d^r^ 



A^s(VP)E ^ e' J ^ J 



27r J {2'kY J (27r)3y (27r)3 7 (27r)3 

xTr [a,SF{E^ + E^, r^, n)a,SF{E2, r^, r^)] . (42) 

The mathematical expression for the subtracted Uehling-diagram can be obtained 
from Eq. ( |^ by the replacement of the bound-electron propagator under the trace 
by the free-electron propagator. 

8 Summary 

In this article, using the Bogoljubov-Parasiuk-Hepp-Zimmermann (BPHZ)-renor- 
malization method, we derived renormalized expressions for the energy shift of 
higher-order for the last two unknown second-order diagrams of bound state QED, 
the two-photon vacuum polarization VPVP and self energy vacuum polarization 
S(VP)E, respectively. It was shown in some detail that counter terms of the outer 
renormalization for the VPVP diagram are cancelled against each other. This result 
simplifies signifiantly numerical evaluations of this diagram. 

In a second part of this paper it has been proven for a wide class of diagrams for 
bound-state-QED which contains the one-loop vacuum polarization VP that the 
counter term of the fourth-rank vacuum polarization tensor vanishes since the sum 
over the quantum number n has been performed until a finite n^ax and as the last 
step in the spherical expansion of the electron propagator. It has been discussed in 
some detail that gauge invariance is not broken in such cases. The disappearance of 
the counter term simplify essentially the evaluation of such diagrams which contain 
the diagram VP. The investigations presented here are to be considered as general- 
ization of corresponding investigations presented in [^ and [|19] . We also mention 
that the numerical evaluation of VPVP and S(VP)E for hydrogenlike ions has to 
be consider as important step in the future to reduce the uncertainty in Lamb shift 
predictions and still is under consideration. Simultaneously, the results of this paper 
establish the use of the BPHZ-approach in bound state QED. 
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Appendix 

The bound-electron propagator is given by 



(43) 



where 9? are solutions of the Dirac equation with Coulomb potential: 

{-ictW + l3 + V{r))^{r)=E^{r) . (44) 

The free-electron propagator of Eq. (^) can be rewritten in a similar form as the 
bound-electron propagator: 



00 2 

S%{E,r,-r,) = jdvY^Y: 



i^s,K,^,{rl)i)s,KA^2) 



(45) 



■^^tl E-Es,p{l-ie) 
where the index s signifies the positive and negative energy states of a free electron 



E,,p = ±^m2 + p2^ (s = l,2) 
and \I/ are the solutions of the free Dirac equation: 

(-zckV + 13) ^(r) = E m{r) 



(46) 



(47) 



The free-electron propagator in spherical expansion can be written as |26 

I G}^^^{E,n,r,)7,\^ Gl%{E,n,r,)nl^ 

- \Gll{E,n,r,)nl' GZ{E,r„r,)rrf 



(48) 



with the radial components (for ri > r2) 



Gl],{E,n,r2) = -{lE + l)^fW+l j\^+y2\-i/2{iri^fW+l) h^. 



(1) 



lK+l/2|-l/2l*^2 



^/W+i) 



Gl^,{E,n,r2) = ~i{E^ + 1) sign(«:) 3\^+y2\^y2{iri^WTl) hfl,/,,_,Jir2^fWTl) 



Gl'],{E,n,r,) = -i{E^ + 1) sign(«:) J|«-i/2|-i/2(^riV^+ 1) hfXy2\-i/2i^'^2VWTl) , 



(49) 
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For the case r2 > ri one simply can employ the symmetry relations: 



In Eq. 



Gl'^{E,r„r,) = Gl],{E,r,,n) , 
Gl],{E,r,,r2) = Gl',{E,r,,n) , 
Gl%iE,n,r,) = Gl'^iE,r,,r,) . 

1^) n^J denote the spin-angualar functions which are given by: 



(50) 



^21 



E xUOi, Vi) xV{d2. ^2) ; vrf = 5: x^«(^i, Vi) X'-I{e2. V2) . (51) 



and Xk is defined as: 



X^(^,0) 



K K,+ ^ — H 



2k+1 






V 



2/. + 1 ^«+ll-iM+|(^'' 



(52) 



/ 



where Y<^_^i,i^^_^i{9,(l)) are the spherical harmonics. The photon propagator reads 
in Feynman gauge: 



D^Luih -t2,ri -ra) 



-g^n 



dE r d^k e"^[-^(*i-*2)-Mr-i-r2) 



2n J (27r)3 E^ - k^ + 



le 



(53) 
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VPVP 



S(VP)E 



Figure 1: The two-photon vacuum polarization (VPVP) and the effective self energy 
S(VP)E. 
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X X 

Figure 2: Potential expansion of the bound-electron propagator. 
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Figure 3: Unrenormalized energy shift due to higher-order contribution in (Za) of 
the VPVP diagram. 
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D E I- 

Figure 4: Lowest-order terms of the potential expansion contributing to the energy 
shift of Fig. I with u = 0. 




Figure 5: Diagram F^ with uj^Ta) = 0. F^ contains a divergent subdiagram 7^ 
one-loop self energy), n denotes a bound electron state. 
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Figure 6: The potential expansion of the unrenormahzed expression Fi^ 
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Figure 7: The inner renormalization according to Eq. (18) 
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aE 



h.o. ren 




Figure 8: This is the final expression for the renormahzation of the VPVP diagram 
if one subtracts the Kallen-Sabry-terms. 
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real photon 

a) b) 




Figure 9: The diagrams of Delbriick scattering a) and energy splitting b). 
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4^ r 3 3^ S 4 4*^ J^ 2 

,G A,H hX, a r h a v_ 






Figure 10: The diagrams of hght-hght scattering. 
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S(VP)E 





Figure 11: The potential expansion of the S(VP)E diagram. The diagrams A3, B2, 
C2 and D are convergent. The diagrams A2, Bl and CI are separately divergent 
but their sum is convergent but not gauge invariant. Diagram Al is the so-called 
Uehling contribution of the S(VP)E. 
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Figure 12: The renormalized expression of energy shift of higher-order of the self 
energy vacuum polarization. 
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